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Abstract. Measurement of redshift-space distortions (RSD) offers an attractive method to directly 
probe the cosmic growth history of density perturbations. A distribution function approach where 
RSD can be written as a sum over density weighted velocity moment correlators has recently been de- 
veloped. In this paper we use results of N-body simulations to investigate the individual contributions 
and convergence of this expansion for dark matter. If the series is expanded as a function of powers 
of /Lt, cosine of the angle between the Fourier mode and line of sight, then there are a finite number 
of terms contributing at each order. We present these terms and investigate their contribution to the 
total as a function of wavcvector k. For the correlation between density and momentum dominates 
on large scales. Higher order corrections, which act as a Finger-of-God (FoG) term, contribute 1% at 
k - 0.015 h Mpc~\ 10% at fc - 0.05 h Mpc"^ at z = 0, while for k > 0.15 h Mpc"^ they dominate 
and make the total negative. These higher order terms are dominated by density-energy density 
correlations which contributes negatively to the power, while the contribution from vorticity part of 
momentum density auto-correlation adds to the total power, but is an order of magnitude lower. For 
fi"^ term the dominant term on large scales is the scalar part of momentum density auto-correlation, 
while higher order terms dominate for fc > 0.15 ft. Mpc^^. For f/' and /i^ we find it has very little 
power for k < 0.15 h Mpc^^, shooting up by 2-3 orders of magnitude between k < 0.15 h Mpc^^ and 
k < 0.4 h Mpc^^. We also compare the expansion to the full 2-d P*^(fc, /^), as well as to the monopole, 
quadrupole, and hexadecapolc integrals of P'*''*(/e, fi). For these statistics an infinite number of terms 
contribute and we find that the expansion achieves percent level accuracy for kfi < 0.15 h Mpc~^ at 
6-th order, but breaks down on smaller scales because the series is no longer perturbative. We explore 
resummation of the terms into FoG kernels, which extend the convergence up to a factor of 2 in scale. 
We find that the FoG kernels are approximately Lorentzian with velocity dispersions around 600km/s 
at z = 0. 
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1 Introduction 

Galaxy redshift surveys are one of the most powerful tools to probe cosmological models [1] . One way 
to extract the information is from the shape of the power spectrum or correlation function, assuming 
it traces the underlying dark matter. Another method involves baryonic acoustic oscillations (BAOs), 
detected in various redshift surveys, which enables one to measure angular diameter distance and 
compare it to the same quantity measured in cosmic microwave background. This in turn probes the 
expansion history of the Universe and allows to study the nature of dark energy [e.g., 2, 3]. Third 
piece of information in redshift surveys comes from redshift space distortions (RSD): the observed 
galaxy distribution is distorted along the line of sight due to the Doppler shifts caused by peculiar 
velocities [4-6]. In linear theory this allows one to measure the rate of growth of structure, which 
allows for another way to measure the matter content of the universe, including the amount and 
nature of dark energy. The last two methods are complementary: cosmological models in different 
gravity theories with the same expansion history cannot be distinguished by the distance scales of 
BAOs, but can if growth of structure is also measured [e.g., 7-9]. Additional information is obtained 
by Alcock-Paczyhski test [10-12]. 

RSD have been analyzed in many galaxy surveys to determine the cosmological models [e.g., 
13-22]. However, it was shown by [23-26] that the growth rate reconstructed from the redshift-space 
distortions can have scale dependent biases, which indicate a breakdown of linear theory predictions. 
These effects show up on relatively large scales, suggesting one must go beyond the linear theory in 
the analysis of RSD. This will become even more important in the future, with several ongoing and 
upcoming galaxy surveys that will measure RSD to a high precision [27-31]. 

Given the high precision of the future surveys, correspondingly more accurate theoretical pre- 
dictions become essential for their interpretation. As was emphasized by [32], there are important 
nonlinear effects that need to be addressed in order to achieve accurate theoretical predictions. In 
order to account for the nonlinearity of the gravitational evolution, standard perturbation theory has 
long been used to describe the power spectrum at quasi- nonlinear scales [e.g., 33]. Recently there have 
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been many studies to predict the power spectrum in nonlinear regime beyond the framework of the 
standard perturbation theory (SPT) [34-39]. These different approaches were compared to the TV- 
body simulations in [40, 41]. Similarly, initial RSD work was based on the lowest order SPT [42-45]. 
However, as pointed out by [32, 43], SPT in redshift space breaks down at larger scales than in real 
space because of nonlinear redshift distortion effects, sometimes called Fingcr-of-God (FoG) effect [4]. 
Recent development using more sophisticated perturbation methods applicable to the redshift-space 
power spectrum includes [46-48]. [22] present detailed comparisons of these predictions with the 
observed galaxy data. 

Recent paper [49] has shown that one can write the Fourier mode of density in redshift space 
as a sum over mass weighted moments of radial velocity, which arc integrals of powers of velocity 
over the momentum part of the phase space distribution function. The corresponding RSD power 
spectrum can be written as a sum over auto and cross-correlators of these moments. This series 
always converges for sufficiently small expansion parameter defined below. We will use the Fourier 
description in this paper and scale is expressed in terms of wavevector amplitude fc, while angular 
dependence is expressed in terms of /i, cosine of the angle between the line of sight and Fourier mode. 
The expansion parameter depends on the product of the two fc|| = kfj,. It has been shown in [49] 
that the moments can be decomposed into helicity eigenstates, which are eigenmodes under rotation 
around direction of k vector. Only equal helicity eigenstates correlate, leading to a specific angular 
structure of the correlators. This analysis shows that if one expands the series into powers of /z^, 
a finite number of terms contribute at each (finite) order. This suggests that RSD can be better 
understood in terms of this expansion rather than the Legendre moments usually used [49]. On the 
other hand, Legendre moments are uncorrelated in real observations, while powers of are not, 
leading to correlations between the higher and lower orders. We will pursue both approaches here. 

This is the second paper in a series studying the redshift-space distortions based on a distribution 
function approach, following the theory and angular decomposition presented in [49]. In this paper 
we test the formalism to describe the redshift-space power spectrum in nonlinear regime using a 
large set of cosmological TV-body simulations, as well as present the individual terms of expansion for 
comparison against each other. We focus on the dark matter in this paper, leaving the application 
to halos and galaxies to future work. The structure of this paper is as follows. In section 2 we 
briefly describe the distribution function approach to RSD. Then we apply it to simulations to test 
this expansion and show the contributions from individual terms: in section 3 we first show the 
contributions to the 2-d power spectrum in redshift space, then proceed to Legendre moments. We 
discuss the FoG modeling in section 4 and present an attempt to compare the expansion to one in 
terms of volume weighted quantities. Finally in section 5 we apply the method to powers of /i^ 
expansion, which we argue is a more natural way to expand 2-d information, showing individual 
contributions to the lowest order terms. This is followed by conclusions in section 6. 

2 Redshift-space distortions from the distribution function 

The exact evolution of coUisionless particles is described by the Vlasov equation [1]. Following the 
discussion by [50], we start from the distribution function of particles /(x, q, t) at phase-space position 
(x, q) in order to derive the perturbative redshift-space distortions. Here x is the comoving position 
and q = p/a is the comoving momentum (p is the proper momentum and a is the scale factor). The 
density field in redshift space is related to moments of distribution function as 




(2.1) 



where H is the Hubble parameter and THk) is the Fourier transform of T.f 



(x), defined as 




(2.2) 



where U|| is the radial comoving velocity, mu\\ = q\\ = q - f, to is the particle mass, f is the unit vector 
pointing along the observer's line of sight and p is the mean mass density. 



- 2 - 



The power spectrum in rcdshift space is given by [49] 

L=0 ■ ^ ^ L=QL'>L ■ ■ \ / 

where k\\/k = cos0 = fj.. It is useful to compare this to Kaiser's hnear theory prediction [5, 32]. Thus 
we have 

f (1 + /A^')'flin(fc) ; linear, 

i^Kaiser(k) " | ^ 2^^^2 p^^ + ^2^4 (^)' Pn 5 uonliuear, ^^'^^ 

where Pun denotes the linear power spectrum and / ~ din D/ din a with D the growth factor. These 
terms will in general have nonlinear corrections, so we call this approximation the nonlinear Kaiser 
order approximation. Replacing these lowest 3 moments with the standard linear theory we obtain 
the original linear Kaiser model of equation (2.4). Here we want to view this series simply as a series 
in fc|| , investigating the convergence as more terms are added. 

Note that the calculations never require anything but simple power spectra of mass-weighted 
powers of velocity to be computed from the simulations. As we will compare RSD power spectrum 
to the sum from individual terms there should not be much sampling variance in the comparison, 
because both are calculated from the same simulation, so the large-scale fluctuations will be the 
same. The order of fc|| = kfj, needed for convergence to a given level of accuracy will inevitably 
increase as one goes to increasingly small scales, with the whole expansion eventually breaking down 
once kfia/H > 1, where cr is a typical (comoving) velocity of the system. We will see in section 3 that 
the nonlinear terms of Pll', particularly P02, have significant contributions even at the scales larger 
than the breakdown scale. 

2.1 Angular dependence 

By performing helicity decomposition [49] show that the power spectrum can be written as 

PLL'{k)= E E (2-5) 

(i=L,L-2,..) {l'=L\L'-2,..; l'>l) m=0 

where P[^{ii) are the associated Legendre polynomials, which determine the angular dependence of 
the spherical harmonics. There are 5 numbers that describe these objects: L and L' describe the power 
of two velocity moments we are correlating, /, I' describe the rank of the object, for example Z = 1 is 
rank-1, which is a 3-d vector, Z = 2 is a 3-d tensor etc. Finally, m is the helicity eigennumber, which 
ranges between and I {I < I'). Only equal helicity components of expansion have a non- vanishing 
correlator. There is a close relation between the order of the moments and their angular dependence. 
The lowest contribution in powers of ^ to P'"'(fc) is ^^"""^ if L + L' is even or '^^ if L -I- L' is 

odd, and the highest is \ Thus for Poo(k) the only angular term is isotropic term (/i"), for 

Poi(k) the only angular term is fj,^, Pii(k) and Po2(k) contain both /i^ and etc. Note that only 
even powers of ^ enter in the final expression, as required by the symmetry. We can thus write 

00 , ^ , s 2L 4L 00 . -.nL /.,xL+L' 2(L+L') 

L=0 ■ ^ ^ j=2L L=QL'>L ' ' \ / j={L+L')oiiL+L' + l) 

(2.6) 

so that terms P[^2' coefficients in expansion in powers of of contributions of L, L' terms to P'^" . 
Since the j index has to be even, the lowest order is either L + L' or L + L' + 1, whichever is even, 
thus the higher order terms also have higher order powers of and one can in principle separate 
them from lower order terms. These terms can be uniquely extracted from simulations from angular 
dependence of Pll' terms and so we will focus on them, although sometimes it is useful to decompose 
them into the individual helicity eigenstates instead. 
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3 Numerical analysis 



3.1 iV-body simulations 

The power spectra of the derivative expansion are all from mass-weighted velocity moments and 
thus can be straightforwardly measured from simulations. We use a series of A'^-body simulations 
of the ACDM cosmology seeded with Gaussian initial conditions, which is an updated version of 
[51]. The primordial density field is generated using the matter transfer function by CMBFAST [52]. 
We adopt the standard ACDM model with the mass density parameter = 0.279, the baryon 
density parameter fit, = 0.0462, the Hubble constant h = 0.7, the spectral index Ug = 0.96, and a 
normalization of the curvature perturvations = 2.21 x 10~^ (at k = 0.02 Mpc""'^) which gives 
the density fluctuation amplitude cs ~ 0.81, which are the best-fit parameters in the WMAP 5-ycar 
data [53]. We employ 1024^ particles of mass 3.0 x 10^^h~^ Mq in a cubic box of side 1600 Mpc. 
The positions and velocities of all the dark matter particles are output at z = 0, 0.509, 0.989, and 
2.070, which are quoted as z = 0, 0.5, 1, and 2 in what follows for simplicity. We use 12 independent 
realizations in order to reduce the statistical scatters. For the detail of the simulations see [51]. 

Next we describe how we measure the power spectra from our simulation samples. We assign 
the density field and the mass-weighted velocity moments in real space on 1024^ grids using a cloud- 
in-cell interpolation method according to the positions of particles. To directly measure P'**(k) we 
also need the density field in redshift space. In measuring the redshift-space density field, we distort 
the positions of particles along the line-of-sight according to their peculiar velocities before we assign 
them to the grid. We regard each direction along the three axes of simulation boxes as the line of 
sight and the statistics are averaged over three projections of all realizations for a total of 36 samples. 
We use a fast Fourier transform to measure the Fourier modes of the density fields in real space 
(5(k) and in redshift space (5s(k), as well as the mass- weighted velocity moment fields in real space 
T|^(k). Finally, the power spectrum in redshift space, P^'^ik), as well as the power spectra of the 
mass- weighted velocity moments Pll'O^), are measured by multiplying the modes of the two fields (or 
squaring in case of auto-correlation) and averaging over the Fourier modes. Throughout this paper, 
we neglect shot noise because we have sufficient number of dark matter particles and such an effect 
is thus negligibly small. Error bars in the following results are estimated from bootstrap resampling. 
The dispersion in power spectra measurements is large on large scales because of sampling variance, 
but it is mostly eliminated by taking the ratio of the two spectra obtained from the same set of 
realizations (e.g. [54]). 

3.2 Matter power spectrum 

In this subsection wc first measure the redshift-space power spectrum, P^''{k,^), directly in redshift 
space, which wc call the "reference" power spectrum. The reference power spectrum in redshift space 
is shown as functions of (fc,/i) at z = and 2 as the red lines in figure 1. We adopt the constant /i 
binning into five bins between < /i < 1, but only three /x bins among the five are plotted. In figure 1 
we also show contributions of the terms of Pll' for (0 < L -I- L' < 4) to P'''^{k, fi) computed from the 
mass-weighted velocity moments. At /i contributions from the higher order power spectra of the 
velocity moments are small and P^'^ ~ Poo because each Pll' is multiplied by a factor of (kfi)'"^^ . 

On large scales one expects Pqq to be followed by the other two linear order terms, which are Pqi 
and the scalar part of Pu, i.e. P^i^. Note however that the latter scales as while there are two 
nonlinear terms, the vector part of Pu, Pn ^, and the scalar part of P02, which itself has two terms, 
both helicity 0, one from energy density correlated with the density Pqq'^, and one from anisotropic 
stress density correlated with density, PqI^i that contain terms which scale as As a result, for 
sufficiently low ji these nonlinear terms dominate over the linear term in PH^ even on very large 
scales. For example, for = 0.1 at z = we sec that P02 dominates over Pu on all scales probed, 
despite the fact that Pu contains a linear order term. We also see that Pu does not follow the linear 
theory on all but the largest scales, but instead has the shape dependence similar to Po2; characteristic 
of the nonlinear terms. As pointed out in [49], the nonlinear helicity 1 (vector or vorticity) terms 
in Pll are closely connected to P02 and partially cancel each other. This angular decomposition is 
discussed further below in section 5, where we present the individual helicity terms separately. 
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Figure 1. Power spectra measured in redshift space P''*(/c,/i) and individual contributions to it from the 
terms of the moments expansion up to 4-th order at z = (left) and z = 2 (right). The width of /i bin is 0.2 
centered around the values shown in each panel. The solid and dashed lines show the positive values, while 
the dotted lines the negative values. 



Because of (k[i)^^^ weight the higher-order terms scale more rapidly with fc, and dominate 
on small scales: this is the region where RSD are dominated by FoG effects. One needs to take 
into account more and more higher-order terms in order to make the expansion (equation (2.3)) 
valid at such smaller scales. This effect is more significant at z = due to higher velocities. One 
can see that the higher order terms cross the lower order terms at k\i ~ 0.2 h Mpc~^ (z — 0) and 
fc/i ^ 0.4 h Mpc~^ (z = 2). This is where the perturbative parameter kjjLa/ H becomes of order unity 
and the perturbative approach breaks down, consistent with [32]. At that point higher order terms 
dominate over the lower order terms and we no longer have a convergence. This can be seen in figure 
1: while for high k the redshift space power spectrum P'"'{k, ji) decreases in power relative to the real 
space case Poo{k), the individual terms in the series expansion increase due to their {kfi)^~^^ . This 
suggests that a non-perturbative approach is needed in this regime: we will explore the so-called FoG 
resummation in section 4. Figure 1 suggests that the typical velocity a in the expansion is about 
500km/s at z = and 250km/s at z = 2. We confirm these numbers in a more detailed FoG analysis 
below. 

3.3 Legendre expansion 

We can compare the agreement between moments expansion and the full P'^'^ as a function of the 
order in the series {L,L'). It is customary to expand the redshift-space power spectrum in terms of 
Legendre multipole moments [e.g., 22, 24, 55-57]. The motivation for this expansion is that if one 
uses full angular information the Legendre moments are uncorrelated. Using Legendre polynomials 
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we have [58] 

p-(fc,/.)= pnmif^)- (3.1) 

1=0,2,4, ■■■ 

The muhipole moments, Pf'', are obtained by inversion of this relation, 

11 + 1 f^^ 

Prik) = J ^ P^%k,ti)Vi{ii)dti . (3.2) 

Most of the studies on modchng redshift-space distortions focus on the monopolc {I = 0) and 
quadrupole {I = 2), ahhough hexadccapole {I = 4) also contains some cosmological information 
[57]. 

In figure 2 we show the monopolc power spectrum at z = 0, 0.5, 1, and 2 summed up to nonlinear 
Kaiser, 2nd, 3rd, 4th, and 6th order approximations in fc/i expansion. Here we denote the summation 
at a given order as including all terms that have the same L-\-L': hence 2nd order includes all 3 Kaiser 
terms PqOj Pqi a-nd Pn, as well as P02, while the nonlinear Kaiser model includes only the first three. 
The lower panels show the error for a given order, Potum/^o rcf ~ 1- ^^'^ linear theory power spectra 
with the input cosmological parameters of our simulations are also plotted for comparison. The power 
spectrum of the nonlinear Kaiser model starts to deviate from the reference spectrum at very large 
scales, k ~ 0.05 h Mpc~^. However, adding the term P02, which has the same order contribution as 
Pii, to the nonlinear Kaiser model, improves the accuracy. Adding the higher order terms continues 
to improve the accuracy down to smaller and smaller scales, but the gains decrease as we approach 
the scale k — (t~^ ~ 0.2 h Mpc~^ {z = 0), where the perturbative expansion breaks down. Our 
formula for the redshift-space monopole spectrum Pg gumj summed up to or more than 6th order, is 
accurate within a few percent accuracy at fc ~ 0.2 h Mpc^^ at z = and at fc ~ 0.4 h Mpc~^ at 
z ~ 2. It predicts not only the overall shape of the redshift-space power spectrum up to these scales 
but also baryon acoustic oscillations (BAO): to see this more clearly we show the summed power 
spectra divided by the smoothed no-wiggle spectrum [59] in figure 3. 

Figure 4 is the same as figure 2, but shows the results for the quadrupole spectra Pj*"*. Because 
the nonlinear quadrupole spectra crosses zero at high-fc, there exists a singularity point for the ratio of 
summed and reference spectra at small scales. The predictions for the quadrupole moment reproduce 
the reference spectrum within a few percents up to the scales of the singular point, k ~ 0.15 h Mpc~^ 
at 2 = and fc ~ 0.3 /i Mpc~^ at z = 2. The quadrupole spectra divided by the corresponding no- 
wiggle spectrum are shown at the right side of Figure 3. Figure 5 shows the results for the hexadecapole 
spectrum. We adopt broader k binning for the hexadecapole moment at fc < 0.1 h Mpc~^ and 
put artificial cuts for the plots of fcP|^ because of large sampling variance. We do not show the 
results obtained from the nonlinear Kaiser, 2nd and 3rd order approximations in lower panels because 
they strongly deviate from the reference power spectrum (as shown in the upper panels of figure 5). 
Although the measurement of the hexadecapole moment of the redshift-space power spectrum is very 
noisy, the higher order expansion predictions give a good agreement if we consider summation to 6th 
order. 

4 Fingers-of-God resummation 

It is clear from the results in previous section that while for kfia/H < 1 we have a convergence 
and only a finite number of terms need to be considered, there is no convergence for k^a/H > 1: 
individual terms become larger and larger as we go to higher orders, yet the total sum in P'^^{k,fi) 
remains well behaved. This suggests we need to explore ways to resum the terms. 

While it is difficult to make progress in general terms, there are specific situations that can be 
controlled. We are interested in a situation where pieces of terms disconnect within the correlation 
function. For example, in P02 term we correlate S with (1 -I- S)up and if is dominated by the small 
scales, as will be the case in systems with large velocity dispersion caused by nonlinear gravitational 
collapse, then on large scales this term becomes P02 ~ PqqO'^, where = and the total con- 
tribution to P'^'^ is —PoQ{kiJ.a/H)^. We see that this term scales as the linear order term, but has 



- 6 - 




Figure 2. Upper panels: we show monopole moments of power spectrum in redshfit space P^°. The vertical 
axis is multiphed by k to clearly illustrate departures from a reference power spectrum. The green, blue, 
yellow, magenta and cyan lines respectively show our model prediction up to nonlinear Kaiser, 2nd, 3rd, 4th 
and 6th order corrections, measured from the simulations. The black lines are linear theory prediction. The 
black points with errorbars show the reference power spectrum. The red lines show our FoG model (section 
4). Lower panels: error between the summed power spectrum and the reference spectrum. The meaning of 
the color of each line is the same as that of the upper panels. For reference errorbars are shown for the result 
of our FoG model. 



opposite sign to it, i.e. This term suppresses power and this suppression scales as fc^. The long range 
correlation is entirely in the density field. While this analysis suggests the mean square velocity field 
enters as the physical parameter, as discussed in [49], any bulk velocity that displaces particles as a 
solid body will not contribute to FoG. This cancellation shows up in part of Pn cancelling P02, such 
that only the dispersion part enters into the term, while the bulk part cancels. 

One can identify similar terms at higher order, for example Po4 ^ ^oo('"|} and one can write 
(wn) = 3(7'' + (un)c, where (wn)c is the connected part of the curtosis. There will also be a term 



P2: 



Pnna . The total contribution to P^'* from these terms is thus 



P^^(fc,^) = Poo[l-(Wff)' + (W^)V2 + 2(W-ff)'(4)c/4!---] 



(4.1) 



The same calculation can be done at the field level. We are interested in the situation where the pieces 
of a term like S{x)uf:{x.) are disconnected within a correlation calculation, e.g., (X{y)5{x)uf:{x 



(X{y)6{x)) (uu)+ other terms. We can re-sum the pieces into a FoG factor 



exp 



(4.2) 



where {u^)c is the connected part of (u^) and note that the odd L terms are zero by symmetry. The 
lowest order term here is just the usual Gaussian kernel exp (— A:^/i^(T^/2iJ ^) . After this re-summation 
we can write the rcdshift-spacc density field as 



Sik) 



y- 

L=l 



ikfi 

IT 



(l + <5(x)) 



(x) 



(4.3) 
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Figure 3. Redshift-space power spectrum divided by the no- wiggle approximation [59], monopole (left) and 
quadrupole (right). The green, blue and cyan lines show our model prediction up to nonlinear Kaiser, 2nd 
and 6th order corrections. The dotted black lines are linear theory prediction, while the solid black lines 
the reference power spectrum. The red lines show our FoG model. The negative values of the quadrupole 
spectrum on small scales are shown as the dotted lines. 



where 



is understood to be minus all possible internal averages of any number of it's. This 
motivates us to write 

P-^^(k) - G([fc/ia/i/]2)PKaiscr(k), (4.4) 

where Pxaiscr account for the lowest 3 terms given by equation (2.4) and where G{x) is exponential 
in the simplest case where higher order reduced moments can be ignored, while more generally it is 
a function with alternating signs of coefficients in Taylor expansion. We have written G{x) in terms 
of a; = [kfj,a/H]^ only: if curtosis is present then we can either write it by adding additional {kji)'^ 
terms to the exponential, exp[— (fctr/i/if)^ + 2(A:/x/7?)^(u||)c/4! • • •], or, equivalently, we allow for a 
more general functional form of G{x) than an exponential. 

Note that in the simplest form this "derivation" gives exactly the exponential FoG form proposed 
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Figure 4. Same as figure 2, but for the quadrupole. The dashed hues at the top panels show positive values 
while the dotted lines show negative values. 
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Figure 5. Same as figure 2, but for the hexadecapoles. We adopt the different bin sizes in logarithmic spacing 
at < 0.1 h Mpc~^ and k > 0.1 h Mpc~^. Artificial cuts are put for the plots of the hexadecapoles at low k 
because of large sampling variance. 



in the literature [12, 32. 60-62]. Other forms for G{x) have been proposed in the literature, e.g. 

Qfj. ^ (f^ /ij~)2\ = / (1 + Lorentzian, 
\ cxp (— x) Gaussian, 



(4.5) 



where a was treated as a free parameter. See [26, 63] for the studies to adopt more than one free 
parameter for the FoG term. These empirical forms all behave qualitively in the same manner, i.e. 
to first order of Taylor expansion in they give G = 1 — {kjia/H)^ + 0{[kiia / H]'^) and higher 

order terms alternate in the signs. 

The effective velocity dispersions can be evaluated from our expression of the redshift-space 
power spectrum (equation (2.3)) and uniquely resummed into the FoG kernels because the higher- 
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order terms of the expansion contain all information. We know that a Taylor series of the FoG kernels 
like in equation (4.5) produce positive and negative terms alternatively, just like the terms in equation 
(2.3). However, there are a lot of Pll' terms in our expansion. At each even order both positive and 
negative terms appear, e.g., at 4th order the Po4 a-nd P22 terms have positive contributions while the 
P13 term has a negative contribution. On the other hand, at odd order contributions are all positive 
and negative in tern, e.g., at 3rd order the terms are all negative and at 5th order all positive. Thus we 
decided to include three independent FoG kernels for the two even linear order order terms (Pqo and 
Pll) and one for the lowest odd order Pqi. These three terms thus need to be individually multiplied 
by something which corresponds to the generalized FoG kernels. To do this, we can apply the same 
reasoning as the discussion around equation (4.1) to other terms that contain long range correlations 
between density and velocity, i.e., portions of terms such as P03 and P13 can be respectively written 
as P03 ^ —QPQi{kiJ.(T/H)^ and P13 ~ —3Poi{kfia/H)^, which, together with the higher order terms, 
can be resummed into a similar function of even orders in powers of kfi/H. 
Finally, we can write 

P^^(k) = Goo {[kfiaoo/H]') Poo + 2Goi {[k^iaoi/H]^) ^Pqi + Gn {[k^^an/H]') ^Pii. (4.6) 

We have defined three different velocity dispersions (Too, <^oij and an. The expressions of the velocity 
dispersions in the FoG kernels can be uniquely derived from equation (2.3) as 

<^oo(k) = ^, (4.7) 

f^oi(k) = 7: p , (4-8) 

^?i(k) = (4.9) 

To go beyond that and determine the form of FoG kernel we expand the FoG terms as a Taylor 
series in terms of (fc/i)^. Because all the phenomenological FoG models have the first derivative equal 
to —1, we define dGLL'{x)/dx\x=o = —1. Now let us consider the ansatz for the FoG model, 

GLL'{xLL';aLL')= (1 + ^) " , (4.10) 

where xll' = [kucTLL' / H)'^ ■ Each FoG kernels contains two parameters, aLL' and aLL'- These 
can reproduce the functional forms in previous studies: Lorcntzian for oll' = 1 and Gaussian for 
aLL' = 00. The FoG parameter aLL' is related to n-th derivative of Gll' as (n > 2) 

^G..,(x)u.oH-irn(i + ^) ■ (4.11) 

m— 1 ^ ^ 

We use the aLL''s determined from the lowest contribution at second order, as 

aLL'{^)=\G'LL'{^)-^\ (4.12) 

Note that the expressions for our FoG models of equation (4.6) preserve full generality up to order 
{kji/H)^ . The quality of the ansatz can be investigated by looking at higher order terms: the com- 
parison of the derivatives of the FoG kernels up to 4th order with the expansion terms is given in 
appendix A. 

The FoG model in equation (4.10) is not the only option to choose. As discussed above, the 
rcsummation can be expressed more elegantly in terms of connected moments, where the cumulant 
theorem naturally leads to an expression of the form 

Gll' = cxp[-(A://(7/i/)2 + [kiiT/Hf + •••], (4.13) 
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Figure 6. Spherically averaged FoG parameters, velocity dispersions (upper panels) and the power law index 
of FoG kernel (lower panels). The circles/triangles have been respectively offset in the negative/positive 
direction for clarity. 

where the lowest two cumulants are variance = and kurtosis = (u|)c, which is the connected 
part of (ujj). We can generahzc this expression and introduce 3 variance and 3 curtosis terms separately 
for 00; 01, and 11 terms, and we can relate the kurtosis terms tll''s to the parameters q;ll''s defined 
above as 2{tll' / (^ll'Y = I/ull'- 

Let us summarize this discussion: the decoupling of small scale velocity dispersion like terms from 
the long range correlations motivates a resummation of the terms into the so called FoG kernels, which 
multiply the long range correlation terms contained in density-density (Pqo); density-momentum (Pqi) 
and momentum-momentum (Pn) correlators. Only a portion of the terms can be motivated in such 
a way, while other terms are simply nonlinear couplings that do not reduce to linear order correlation 
on large scales. It is therefore difficult to provide a formal justification for this resummation, but 
it is worth analyzing to what extent this approach is useful. Here we use the distribution function 
expansion in equation (2.3) to formally define FoG kernel parameters: up to (kfi)'^ this is equivalent 
to the exact expansion, but requires 3 different FoG kernels. An additional point that needs to be 
emphasized is that these quantities as we defined are a function of angle and scale, i.e. they are fi 
and k dependent. Below we use numerical simulations to compare the original expansion in equation 
(2.3) to the FoG resummation version to see to what extent FoG approach is useful for the general 
treatment of RSD. 

4.1 Testing the Fingers-of-God model 

In this subsection we compare the FoG resummation of the higher order terms discussed above to 
numerical simulations. The upper panels of figure 6 show velocity dispersions determined from our 
simulations using equations (4.7) - (4.9). For clarity we plot the spherically-averaged velocity disper- 
sions, i.e., those obtained from the monopole spectra. While the 3 velocity dispersion terms are similar 
to each other, they differ in the amplitude and scale dependence, suggesting that it is important to 
independently consider the FoG kernel for each spectrum in the nonlinear Kaiser formula. Note that 
these values are higher than those determined by previous studies (e.g., figure 7 of [47]). This is 
because our quantities are based on the mass-weighted velocities in contrast to the volume-weighted 
velocities discussed in [47]. We discuss the latter in section 4.2. 
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Figure 7. Angular dependence of FoG parameters at z = 0. The width of \x bin is 0.2 centered around the 
values shown in the top left panel. 



The lower panels in figure 6 show the inverses of the seeond parameters for our FoG model, 
l/o^ii'l^) defined by equation (4.12). As we have discussed above, the parameters are equivalent 
to 2 {tlli / ffLLi) when we adopt the kurtosis as the second parameters. Note that for a Lorentzian 
model and a Gaussian model we have I/qll' = 1 and l/aLu = 0, respectively. The impression from 
the figure is that a parameters have strong k dependences and behave differently from each other at 
large scales. However, the FoG terms have negligibly small effects on the shape of the power spectrum 
at large scales. On the other hand, at lower redshifts the a values converge to nearly a constant at 
k > 0.1 h Mpc~^ where the FoG effect starts to play an important role, with typical values between 
1 < otLL' < 2, i.e. Lorentzian FoG kernel is a better approximation than Gaussian. The convergence 
of a^L' to a single value is worse as we go to higher redshifts and the difference between the values of 
dLL'^^ remains large. In order to see the angular dependence of our FoG parameters, we show ctll' 
and 1/aLL' at z = as functions of k and at the upper and lower panels in figure 7, respectively. 

Now we focus on the power spectrum P*'* using the FoG model discussed above. In figures 2. 4 
and 5 the resulting power spectra with the FoG kernels are shown as the red points. As we have seen 
in section 3.2. there were certain scales at which the power spectrum from our perturbative expansion 
approach breaks down and diverse, k ~ 0.2 h Mpc^^ at z = and k ~ 0.4 h Mpc~^ at z = 2 for the 
monopole spectrum even though we sum up to the 6th order terms. Our FoG model dramatically 
improves the results. The accuracy of a few percent is achieved up to fc ~ 0.4 h Mpc~^ at z = and 
even up to to fc < 1.0 /i Mpc~^ at z = 0.5 and 1. On the other hand, the improvement at z = 2 is not 
so much better than the lower redshfits because of a strong suppression by the FoG. Our FoG model 
works for the quadrupole spectrum as well as for the monopole spectrum: it predicts the quadrupole 
spectrum down to very small scales (see the upper panels of figure 4 and the right panels of figure 3). 
Even the accuracy of the hexadecapole spectrum of our FoG model is quite good, at the level of 10%, 
and it reproduces the shape over all the scales probed (see upper panel of figure 5). 

4.2 Mass- weighted vs volume- weighted velocities 

It is worth comparing our approach based on the power spectra of mass-weighted velocity moments 
with those of volume- weighted velocity moments, which are commonly used to model nonlinear power 
spectra (e.g., [25, 32, 47, 62, 63]). The difference comes from the fact that the L-th moment of the 
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mass-weighted moments T^^ = (1 + d)u^ in equation (2.1) contains contributions not only from L-th 
order in perturbation theory, but also {L + l)-th order because of the term (S(x)u|'(x). By regrouping 
equation (2.1) with the same order term, we obtain 



. 7 OO 

6s{k,n) = 6{k,^i) + — U||(fc,pi) + 



ik\\ 



(L-1)! V n 



L-1 , / 1 \ L 



1 f 



k LWH 



(4.14) 

where the first and the second terms of the right-hand side are respectively the zeroth and first order 
terms, while the bracketed terms correspond to the L-th order terms (L > 2). By squaring equation 
(4.14), we obtain the same equation as equation (2.3), but the terms regrouped into the same order 
in SPT, 



P-(k)=Poo + 2^Po«S+%^^, 



Uu u 



+ ••• (4.15) 



where w^^ = — u^^ = [fejj^jk, thus P^.^l = Pol — ^ou^ similarly for higher orders. 

As many previous studies have already discussed (e.g. [63]), measuring the power spectrum of 
volume-weighted velocity moments is not as straight-forward as measuring mass-weighted velocity 
moments used in our formalism. In order to measure the moments of volume-weighted velocities, we 
divide the interpolated moments of mass-weighted velocities by the interpolated density before the 
field is Fourier-transformed T'||'/(l -I- (5) = u^y This can be noisy: some points on the grid may not 
have any particles, so a sufficiently coarse grid is needed. More importantly, the results depend on 
the grid size, i.e. on smoothing (see [32, 64] for a detailed discussion of how to measure the volume- 
weighted velocities from A^-body simulations). Following the same process as described in section 
3.1, wc compute the redshift-space power spectrum based on the power spectra of volume- weighted 
velocity moments up to a given level of accuracy. 

The top set of figure 8 shows the ratio of the summed power spectra of equation (4.15) to the 
corresponding reference spectra in redshift space. Because equation (4.15) is essentially the same 
as equation (2.3), here we want to see the convergence of these expressions. At the order of first 
3 (Kaiser) terms, the expansion with volume-weighted velocities is somewhat closer to the reference 
spectrum than that with mass-weighted velocities, but both approximations are bad. The bottom 
set of figure 8 shows the same quantities as in the top set but those at fc = 0.117 h Mpc~^ at from 
3rd to 6th orders. At higher orders the convergence is slightly faster with mass- weighted quantities 
for monopolcs when they are compared at the same order: one can predict well the redshift-space 
power spectrum at /c = 0.1 h Mpc~^ by including the 4th and 3rd order corrections at z = and 
z > 0, respectively, for our power spectrum of mass-weighted velocity moments, while one needs 
to include the 4th order corrections for the power spectrum of volume- weighted velocity moments. 
For the quadrupole spectra, the situation is more complicated but the value of the spectra from the 
volume- weighted velocity moments expansion is deviated from the reference value at 3rd order at very 
large scales as seen at the bottom panels of the top set in figure 8. Both approaches break down once 
one enters the non-perturbative regime (fc > 0.2 h Mpc~^ at z = 0). We conclude that in terms of 
rate of convergence, there is no advantage in defining volume weighted quantities. If one works with 
galaxies and halos number density weighting becomes essential, since it is difficult to define volume 
weighted velocity moments in a sparsely sampled system [63]. In that situation our approach is the 
most meaningful way to define physical quantities that enter in RSD description. 



5 Expansion in powers of fj? 

As we have seen in previous two sections the series expansion of equation (2.3) is convergent on large 
scales, but not on small scales. For sufficiently high fc any finite order summation fails drastically. 
FoG resummation approach fares better, but even that fails for high fc. One can sidestep these issues 
by considering an alternative expansion in terms of powers of /i^: as discussed in section 2.1 and in 



- 13 - 



7 0.04 

0) 

CO 

w o 

9; 

E 

Qf-0.04 
7 0.04 

Q; 

E 



rj' 

. MonopoleJ- ^ 

NL Kailser 

■ 2nd \ 

3rd , \ 


iii 

/ z=o.5 : 
/ 

■ 4th \ 

6th \ 

1 1 J. 1 1 1 1 1 1 1 Mill 


— - 

\ \ 
\ 

1 — 1 1 ^ 1 1 1 1| 1^1 INI 


/ ' ' " 

/ 7 
z=2 /" 

\ - 


1 1 1 1 1 1 1 1 H 1 1 1 1 1 

■ Quadrupole\ 
1 

" 1 I : 


\ \ 
\ \ 

- : 


r ^ " 

1 


— 1 — 1 1 1 1 1 1 1|. — 1 — 1 1 ii 1 1 

■ .^^-^ -\ " 
y V ' 

. 1' : 

■ ^ W : 
1 . .\j . . . 



0.01 0.1 0.01 0.1 0.01 0.1 0.01 0.1 

k[hMpc"''] k[hMpc"''] k[hMpc"^] k[hMpc"''] 




34563456345634563456345634563456 



Pq P2 Pg P2 Pq P2 Pq ^2 



Figure 8. ( Top set) Error between the summed power spectrum in terms of volume weighted velocity moment 
expansion and the reference spectrum. Lines are the same as the lower panels of figures 2 and 4. (Bottom 
set) Error between the summed power spectrum and the reference spectrum at = 0.117 h Mpc~^. The 
left and right sides of each panel show the monopole and quadrupole spectra, respectively. The horizontal 
axis is a given order in terms of mass-weighted moment expansion (red points) and volume-weighted moment 
expansion (black points). 



[49] for any finite power of /z^ there is a finite number of Pll' terms contributing to it. For /i° only 
Poo contributes, for /x^ Poi, Pn and P02, for ^'^ Pn, P02, P03, P12, Po4, P13 and P22 etc. 

Only these 3 lowest terms, /i", /x^, and /i**, contain cosmological information at the linear order, 
so in principle these are the only relevant terms. However, if we expand the full P*'*(fc, fj,) into powers 
of /x^ and try to determine the coefficients from the data, the resulting coefficients will be correlated: 
only Legendre expansion assures uncorrelated values. As a result there will be mixing of higher powers 
of fi^ into lower powers if they are not accounted for in the fits, or there will be strong degeneracies 
and the fits will be unstable if all the coefficients are accounted for but we allow them to take any 
value. Typically one solves this by regularizing the expansion, i.e. by constraining them to a certain 
range of values. In this paper we will not focus on methods how to determine the coefficients of such 
expansion from the data, but we will show fi^ and /x^ expansion terms to develop some understanding 
of their scale dependence and amplitude. 

5.1 /i^ terms 

In the top panels of figure 9 we show these individual term contributions to the lowest order powers 
of /i (we do not show /i'^ term, which is just the usual real space power spectrum Poo). For /i^ we see 
that the Poi dominates for low k, as expected since that is the only term that does not vanish in linear 
theory. This term follows linear theory prediction for low k, while for k > 0.1 h Mpc^^ it exceeds it, 
just like it happens for the dark matter power spectrum Poo itself. This is not surprising: as shown 
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Figure 9. Contributions to (upper panels), /i'' (middle panels), and /x® and terms (lower panels) in 
redshift-space power spectrum P'"'{k,fi). The positive contributions are shown as the solid and long-dashed 
lines while the negative contributions as the dotted lines (linear theory predictions are shown as the short- 
dashed black line). 

in [49] we can write Pgf (//^) = '^^""/q ; so this term is given by the time derivative of the dark matter 
power spectrum and has a similar scale dependence relative to the linear power spectrum. 

The next term in terms of relevance is Po2- This term has contributions from the correlation 
between the energy density and the density Pqq'^, as well as from the scalar part of the anisotropic 
stress correlated with the density Ppif, both helicity scalars. As discussed in section 4, we expect 
the first term to be dominant and scale as —k^Pooa^/H^, hence to dominate over Pqi for ka/H > 1. 
This gives a ~ 500km/s at z = 0, decreasing to 200km/s for z — 2. As expected we see this term is 
always negative. 

The third term contributing to fi^ is the vector part (helicity 1) of auto-correlation of momentum 
density with itself Pll^. This term is always positive and partially cancels Po2- As discussed in [49], 
it cancels all the bulk motion contributions to cr^. While this term scales in a similar way as P02, it 
is 3-10 times lower in amplitude, so it cancels only a small part of P02. 

We see that the total sum never exceeds the linear power spectrum and becomes negative for 
k ^ 0.17 h Mpc~^ at = and k 0.5 h Mpc~^ at z = 2. For scales smaller than that the /i^ term 
is negative as a consequence of a strong FoG effect. 

5.2 /x^ terms 

This term receives contributions from 7 different terms, Pn, P02, P03, P12, Po4, P13 and P22. They 
are shown at the middle panels in figure 9. On large scales the dominant term is Pn which contains 



- 15 - 



a linear order contribution. This term agrees with hnear theory prediction for k < 0.1 h Mpc~ and 
is above that for k > 0.1 h Mpc""'^, just hke in the case of Pqo and Pqi- 

The next order term in significance should be Po2- We see this term is relatively small and 
does not dominate anywhere. P02 contribution to /i^ arises entirely from from the scalar part of the 
anisotropic stress correlated with the density Pq2'^, which contains /i*^ term contribution to P'*''*, while 
Pqo" docs not. Physically one expects the small scale velocity dispersion to be relatively isotropic, 
hence anisotropic stress density should be small compared to the energy density. 

Next two terms in terms of L + P' arc P03 and P12. These two terms provide the dominant 
correction to Pn on intermediate and large scales k < 0.1 h Mpc~^. On very small scales terms P04, 
P22 and P13 dominate. As discussed in FoG section we expect the first two terms to scale as 3fc^cr''Poo 
and be equal in amplitude, while P13 should be negative and cancel out the bulk flow part of the 
other two terms. We see that this expectation is borne out in simulations: in total these terms add 
power on small scales. The transition happens at a similar scale as for /i^ term, k ^ 0.17 h Mpc~^ at 
z = and k ^ 0.5 h Mpc"^ at z = 2. 

5.3 Higher order terms 

At order higher than jj,'^ we do not have any linear order contributions, so we expect these terms to be 
small on large scales. There are many terms that contribute. At the bottom panels in figure 9 we show 
the total contributions to /.i® and /i® terms. We can see that these terms are indeed negligibly small 
ai k < 0.1 h Mpc~^ at z = and k < 0.3 h Mpc~^ at z = 2. At smaller scales these contributions 
increase with the scale dependences of k^ and k^, respectively. 

6 Conclusions 

In this paper we used numerical simulations to investigate the distribution function expansion ap- 
proach to the redshift space distortion power spectrum [49] . The power spectrum in redshift space can 
be written as a sum over correlators between mass-weighted velocity moments. We analyzed a large 
set of cosmological A^-body simulations to test how accurately this formalism predicts the true power 
spectrum in redshift space and how many terms are needed to achieve a given precision. We compared 
to the RSD power spectrum as a function of wavevector k and cosine of the angle between the Fourier 
mode and line of sight /i, as well as to the lowest Legendre moments, monopole, quadrupole, and 
hexadecapole. These comparisons revealed that the expansion is accurate within a few percent up to 
k ~ 0.2 h Mpc~^ at z = and k ~ 0.3 h Mpc~^ at z = 2, if the corrections up to the 6th order are 
taken into account. All expansions break down at higher k, where the expansion parameter becomes 
larger than unity. 

The expansion can be resummed into power suppression factors known as the Fingcr-of-God 
(FoG) kernels. Our expansion formula suggests that one needs three independent FoG kernels, each 
of which is multiplying the lowest 3 terms in the expansion, which are the 3 terms that contain linear 
order contributions (Kaiser formula). We have found that the 3 velocity dispersions differ in their 
values. We also found that the shape of the FoG kernels differ and are a function of scale, redshift 
and angle. There is no single FoG kernel that would fit the simulation data. This FoG model has 
validity comparable or better than the highest (6th) order summations we have tried. Our FoG model 
predicts the monopole power spectrum with a few percent accuracy up to fc ~ 0.4 h Mpc^^ at z = 
and is accurate for fc < 1 /i Mpc~^ at z = 0.5 and 1. 

The difficulty of obtaining convergence or a perfect FoG model on small scales in Legendre 
moments expansion do not exist in the alternative expansion in powers of /i^-' . This expansion has 
the advantage of containing a finite number of terms for a given value of j, and hence no infinite sums 
are needed, there is always convergence for any finite j. In contrast every Legendre moment receives 
infinite number of terms. We have used the angular decomposition of [49] to derive the individual 
terms for the lowest values of j. We compared their amplitudes as a function of scale, showing that 
for and fi^ the nonlinear FoG terms dominate the Kaiser terms for k > 0.2 h Mpc~^ (z ~ 0). 
Similarly, higher orders (fj,^ etc), which do not receive any contributions at the linear order, are nearly 
zero for k 0.2 h Mpc~^, but grow very rapidly on smaller scales. Because the higher order terms 
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of Pll' also have higher order powers of /i^, one can in principle separate them from lower order 
terms, put them all into terms that we do not care about, and we can marginalize over. Note that the 
expansion in powers of has the disadvantage that the coefficients extracted from observations will 
be correlated, unlike Legendre moments, which are orthogonal, at least on scales small compared to 
the size of the survey. In practice this means that some priors will have to be assumed for the terms 
with higher powers of fp, at least at high k, otherwise they will contaminate the determination of /i*', 
fi^ and /i^ terms. In this sense we cannot escape the high order FoG complications. However, it is 
not clear there is any useful information remaining at fc > 0.2 h Mpc^^ (z = 0) due to the nonlinear 
evolution anyways. This will be explored further with perturbation theory in paper III of this series 



We also compared our approach based on the power spectrum of mass- weighted velocity moments 
to the one based on the power spectrum of volume-weighted velocity moments. We found that the 
power spectrum based on the mass- weighted velocity moments converges slightly faster, although the 
difference is not so significant. The fact that there is no natural way to expand RSD into volume 
weighted quantities, and that using volume weighted quantities does not improve the convergence, 
implies that there is no good reason to work with volume weighted quantities in RSD studies. 

This advantage of mass weighting quantities becomes crucial when the same formalism is applied 
to galaxies and halos, where volume weighted moments cannot be easily defined. It is difficult or 
impossible to measure the velocity power spectrum from the sparse density field and the results 
depend on the smoothing (see e.g. [63] for an attempt and failure to measure it). In contrast to 
this, the analysis presented in this paper can be naturally extended to the clustering analysis of halos 
and galaxies, since density weighted moments are not affected by sparseness of the sample. Recently 
it was shown by [24] (but see also [66]) that the RSD parameters reconstracted from redshift-space 
distortions of dark matter halos have strong halo-mass and scale dependence even on large scales 
(fc < 0.1 h Mpc"^). This can be a serious problem when one wants to constrain dark energy or 
modified gravity theories. This complicated scale dependence is likely to be a consequence of the 
number density weighting of velocity moments, which differs from the mass weighting of dark matter 
if galaxies or halos have a different spatial distribution. The scale dependence will arise even for a 
linear bias model, except on very large scales, where galaxy overdensity Sg <^ 1 [49, 67]. The formalism 
presented in this paper can be used to investigate this scale dependence and we plan to present this 
analysis in paper IV of this series. 
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A Derivatives of FoG kernels from higher-order Pll''s 

In this appendix, we present derivatives of the FoG kernels: Gqo, Gqi and Gn for the auto power 
spectrum of the redshift-space density field, 
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Note that these equations are general and not dependent of the ansatz of the kernel (equation (4.10)). 
As discussed in section 4, the first derivatives of Gll' give the definitions of the velocity dispersions 
(Tll' • The power low index parameters of our FoG model, a^L' , are obtained by substituting equations 
(A.4) - (A.6) into equation (4.12). 
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